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SUMMARY 


A local linearization technique has been successfully used to solve a 
system of stiff differential equations associated vd.th a Magnetic Bearing 
Assembly. The technique has proven to be accurate, stable and extremely 
efficient. A variable order Adams method with a stiff option was used as a 
reference case. 


INTRODUCTION 

Certain models encountered in simulation studies require the solution of 
systems of stiff, nonlinear differential equations. Solving a system for the 
state as a function of time using general purpose flexible integration 
routines becomes costly because of overhead, satisfying error tolerance 
requirements and the fine integration step size which may be required to 
represent properly the system response. The problem is compounded when one or 
several of these models are required to be interfaced as subsystems of a large 
dynamic simulation which may have many time consuming derivative evaluations 
as well as slowly varying states. In such applications, the stiffness of the 
entire system may be amplified, thus requiring the entire simulation to 
operate at a very small integration step size. A method for solving this 
problem is to apply a stiff differential integration technique to only the 
stiff models of the simulation, thus allowing for stiff and nonstiff 
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integration algorithms to be used, and consequently improving the efficiency 
of the overall simulation. Therefore, a numerical integration technique, 
preferably with one derivative evaluation per step, which is stable, accurate, 
and efficient is needed to solve these stiff simulation models. 

An integration technique which exhibits the desired characteristics is 
the local linearization (LL) technique described in reference 1. The LL 
algorithm is especially suited for models for which the Jacobian matrix is 
easily obtained. 

An excellent model to demonstrate the application of the LL technique is 
the Magnetic Bearing Assembly (MBA) described in appendix A. The MBA is a 
magnetic actuator which is used in critical solar, stellar and earth pointing 
applications on Space Shuttle missions. This model has eigenvalues ranging 
from 10 to 2000. However, if the MBA's are interfaced into the Annular 
Suspension and Pointing System (ASPS) simulation model, eigenvalues ranging 
from 0.1 to 2000 are obtained, thus considerably amplifying the stiffness of 
the system. The ASPS is described in reference 2 

The purpose of this study is to demonstrate the LL technique to simulate 
a MBA or other similar models which cannot be solved efficiently using general 
purpose, flexible integration methods. 

The LL performance is documented in this report against a general 
purpose, variable order Adams method with a stiff option (VOADAM) which 
assumes a dense Jacobian. 


B, C, u, V, r, q, 6 
^1 ’ ^ ’ ^3 ’ ^4 


,\d» 



SYMBOLS 

submatrices required for matrix Inversion 
MBA current loop parameters 
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'^AC’ 

A(t) 

"ij 

Cp C2, C3 
D 

? 


f 


U 


f 


L 


f 


f 


c 


F 


max 


h 

I 

ic 

lo 

K 


m 

P, Q 


s 
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MBA circuit parameters 

time varying Jacobian matrix 

elements of B, where i, j = 1, 2, 3, 4, 5 

elements of A(t), where i, j=l, 2, 3, 4, 5 

MBA scale factors 

matrix representing special 5x5 inverse 

elements of D, where 1, j = 1, 2, 3, 4, 5 

n-dimensional vector composed of general nonlinear 

time-varying functions of the state vector Z(t) 

force generated by upper magnetic pole, N 

force generated by lower magnetic pole, N 

MBA force output, N 

force command input to MBA, N 

maximimi MBA force output, N 

actual magnetic gap, m 

nominal magnetic gap, m 

integration step size, seconds 

identity matrix 

MBA current command, A 

MBA bias current, A 

•y 

MBA force constant, Mn/A^ 
payload mass. Kg 

matrix coefficients in general LL algorithm for 
solution of nonlinear time-varying systems 
Laplace operator 
time, seconds 


rri rri frt 

if, J-2’ ^3» ^4» ^5 


common terms used in inverting the 5x5 matrix B 



V 

z(t) 

A 

^8 

li, TI3, 

“b 

C 

A 


MBA voltage conmand, V 

n-dimensional vector representing system states 
difference between nominal and actual magnetic gap, m 
estimate of Ag,m 

functions of magnetic gap and gap coefficients 
MBA bandwith, rad/ sec 
actual gap coefficient, /m 
estimated gap coefficient, /m 


Subscripts : 
k 
L 
n 


integer constant 
denotes lower MBA pole 

denotes steps in construction of submatrices, n = 2, 


3 , 4 , 5 

U denotes upper MBA pole 

Dot over symbol denotes time derivative 
Abbreviation: 


ASPS 

CPU 

LL 

MBA 

VOADAM 


Annular Suspension and Pointing System 

Central Processing Unit 

Local Linearization 

Magnetic Bearing Assembly 

Variable Order Adams with stiff option 
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PROBLEM DESCRIPTION 


Shown in figure 1 is a block diagram of the MBA. The simplified ASPS 
dynamics shown was used purely for MBA model checkout and is not included in 
the state equations. The state differential equation of the MBA is given by: 

Z = F (Z, f^, Ag, Ag, Ag, t) (1) 

A 

where Z, f^, Ag, Ag, Ag and t represent the state vector, the command 
force, gap, gap rate, gap estimate, and time, respectively. The LL solution 
at is given by: 

\+l " \ \ \ \ FF ^c^’ ^®k’ ^®k’ ^®k’ *^k^ 


where 


P 




I) 


(3) 


Q = a“^ (P - hi) 


(4) 


A = 3Z/3Z 


(5) 


— — — 

3F _ 3F c , 3 F 3Ag 3 F 3Ag 3 F 3 Ag 

3t 3f 3t 3Ag 3t ... 3t / 3t 

c ° 3Ag 8Ag 


( 6 ) 
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I = nxn identity matrix 


and h is the integration step size. The order of the system is n, which for 
this system is 5 for the upper pole and 5 for the lower pole of the MBA. The 
n-dif ferential equations of (1), the A matrix (5) and the n-vector (6) are 
given in appendix B. 

The approach taken to efficiently solve equation (2) was to substitute a 
first order Fade " approximation for the matrix exponential and to develop a 
special matrix inverse to take advantage of the zero elements of the resulting 
5x5 matrix. Another approach which was not investigated, but should also be 
efficient when used in conjunction with the first order Fade " approximation 
is to solve systems of differential equations in lieu of computing inverses. 
This would probably reduce the number of operations required to obtain a 
solution as compared to the approach taken. 

If the Fade ^ approximation 

Ah —1 

e == (I - Ah/2) (I + Ah/2) (7) 

is substituted into equations (3) and (4) F becomes: 

F = h (I - Ah/2)“^ (8) 


and Q becomes : 


Q = Fh/2 


( 9 ) 


The matrix 



can be inverted using standard system matrix 
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routines. However, to take advantage of the many zero elements, a special 
inverse was obtained using the partioning and bordering technique discussed in 
reference 3 and summarized in appendix C. This final code resulted in 
approximately a 76 percent reduction in execution time required for the 
inverse problem. 


RESULTS 

Several test cases were examined comparing the LL and VOADAM results. 

Forcing functions used were a step, triangular wave and sine wave. The 

triangular wave had a frequency of 1 Hz and a slope of 1 Newton/second. The 

sine wave had an amplitude of 4 Newtons and frequencies of 2.5 and 10 Hz. The 

LL test cases were run with a integration step size of 0.001 seconds. The 

error criteria used for the VOADAM solutions was selected such that the 

solution was to have no more than 125 units of error per 1,000,000 units of 

magnitude. Two test cases were made with each forcing function except for the 

step input where only one case was run. The two cases run were with and 
3F 

without the -t— vector zeroed. The step, which was input at t = 0.005 seconds 

O L 

3F 

with a magnitude of 1.5 newtons, was only run with zeroed since by 

O t 

definition the derivative of a step has infinite slope. The purpose of 

3F^ 

examining the effect of zeroing the r— ■ vector is that in practice it is 

O L 

seldom available. In some applications a backward difference calculation of 

3 ¥ 

can be used, however, caution should be exercised when using this 
procedure. 

The LL and VOADAM solutions to a 2.5 Hz sine wave is shown in figure 2. 


The variables plotted are the MBA output, f, the error between the VOADAM and 
LL solutions, and the magnetic gap change, Ag, resulting from applying f 
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to a 90 Kg mass. Figures 3 to 5 show the LL and VOADAM solutions to a 10 Hz 
sine wave, a triangular wave and a step input, respectively. Presented in 
tables I to IV is a digital representation of portions of the data shown in 
the time history plots. The maximum error which occurred during the run is 
shown enclosed in the box. 

The data presented shows LL to be extremely accurate even with the 
3 F 

vector zeroed. As mentioned previously, the LL data presented was run 
0 t 

at only one Integration step size, 0.001 seconds. No attempt was made at the 

writing of this report to vary the step size or to investigate methods of 

8 F 

approximating the t — vector with backward differences or other numerical 

o t 

methods of computing derivatives. It is reasonable to assume, however, that 

3 F 

the LL results for — equal to zero could be improved by one or both of the 

O t 

above methods. 

Computer execution times comparing LL to VOADAM for an input sine wave 
frequency of 1, 5, and 10 Hz are shown in the following table: 



1 

5 

10_ 

LL 

2.7 

2.7 

■ 

VOADAM 

24.2 

31.1 

28.0 

%CPU saved 

88.8 

91.3 

90.3 


The execution times are expressed as CPU-second/second of run time and the 
percent CPU saved is the amount of CPU time saved by LL. For timing purposes, 
a step size of 0.001 seconds was used for LL while VOADAM adjusts the step 
size to satisfy a maximum error criteria of 125 units of error per 1,000,000 
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units of magnitude. VOADAM was required to return to the calling program at 
0.01 second intervals. Although this is not the most efficient use of VOADAM, 
it was required to simulate the operating environment of the MBA. The MBA 
receives inputs at 0.01 second intervals from a digital controller. 


CONCLUDING REMARKS 

A local linearization technique (LL) has been successfully used to solve 
a system of stiff differential equations associated with a Magnetic Bearing 
Assembly. The LL technique proved to be accurate, stable and extremely 
efficient when compared against a general purpose flexible Adams integration 
method with a stiff option. In large dynamic simulations, which require the 
simulation of models such as the Magnetic Bearing Assembly, the LL technique, 
when applicable, appears to be one of the most efficient methods to employ. 
When considering budget constraints, the LL technique may very well determine 
whether or not models such as the MBA can be included, without simplication, ’ 
in large dynamic simulations. 

The LL technique presented in this report is not restricted to the 
application discussed herein. For example, in the linear time- varying case, 
the Jacobian and LL coefficients (P and Q) are constant, which makes the LL 
technique especially suitable for real time simulation applications since only 
one initial computation of the Jacobian is required. 
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APPENDIX A 


DESCRIPTION OF THE MAGNETIC BEARING ASSEMBLY 
The Magnetic Bearing Assembly (MBA) Is a magnetic actuator used for fine 
pointing control in the Annular Suspension and Pointing System (ASPS) which is 
described in reference 2. Shown in figure 1 is a block diagram of the 
hardware model and the associated current loop electronics. Actuators of this 
type are inherently nonlinear. This particular model uses bias current 
linearization to remove the current squared nonlinearity and a signal 
proportional to the gap is used to multiply the coil currents to compensate 
for the inverse-gap-squared relationship. A current control loop is used to 
obtain the desired actuator response. The MBA parameters are shown in table 
V. The current control loop parameters are tabulated in table VI. The 
following functional relationships are used to calculate the bias current and 
current control loop parameters. 


I 

o 



(A-1) 


\d ’^LE ^g “b 


(A-2) 


K 

P 


^^®o ^AC ^®o S ''' ^LE ^AC^ “b 


(A-3) 


K 


I 


^AC^ 



(A-4) 



(A-5) 
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(A-6) 


where is the bandwidth and is the maximum force output. The 

parameters Kj^ and K2 were selected to give desired damping 
characteristics. The value of the remaining parameters are the result of 
laboratory measurements. 
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APPENDIX B 


STATE DIFFERENTIAL EQUATIONS, "A" MATRICES AND THE VECTORS 

a t 

The differential equations are 


Upper Pole 


^lU “ h ^^CU “ ^4U^ 


(B-1) 


^2U ^^lU ^1 ^^CU ^4U^ ~ ^4 ^2U^^^3 


(B-2) 



(B-3) 


^4U L ^5U \c “ ^4U \c “ ^ ^4U^ 

Lr« 


(B-4) 


• ho C„ Ag 

^5U " Ii“ \c ^4U 2 " ’^AC ^5U^ 

o ho 


(B-5) 


Lower Pole 


^IL = ^ (^CL -^4L^ 


(B-6) 
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(B-7) 


where 


Z2L = - "4L> - ^4 


"3L = "2L 


(B-8) 


Z,. = (V^ + R.^ - Z,, - IL Z,,) 

4L L^„ L 5L AC 4L AC D 4L 

LE 


(B-9) 


• C-Ag Zc, 

ZcT = T— (+ R»o Z,- + X R. _ Z„ ) 

5L Lg AC 4L 2 AC 5L 

■ o h/ 


(B-10) 


^CU = "l> 


(B-11) 


^CL = S> ^2 


(B-12) 


Tij^ = (1 - CAg) 


(B-13) 


ri2 = (1 + ^Ag) 


(B-14) 


= (1 - cAg) 


(B-15) 
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= (1 + cAg) 


(B-16) 


^2U ^2U ’*' ^3U “'^LIM - \ - \lM 

\ - >4jD ^2L ^ ^ ^2L ^ "3L "''lIM < \ < ''lIM 


s ■ h 


^2 ~ 


F 


U 


(€3/113)2 


2 

5U 


^L - 4 



The 8F/3t vector for the upper pole is given by: 


(B-17) 


(B-18) 


(B-19) 


(B-20) 


(B-21) 


(B-22) 


(B-23) 
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3Z 


lU 


9t 


= K2 [n 


. 3fc 

1 ^1 at 


- (I + c, 




at 


] = K, 


VARl 


^ - iq- 



0 


'^leS 


. VARl 


az 


4U 


= 0 -V > V > V 

at LIM ^ U LIM 


az 


5U 1 


at Lg^ 


^AC^^5U “ 


S^f_f5u_A Ms. ^2 ^5U aAg 


at 


TI3 at 


The aF/at vector for the lower pole is given by: 


az 


-^= K„ 

at ^ 


-C n — +(I - C f)c — ® 

^2 at ^ o ^1 ^ at 




. VAR2 


^^2L h 
— = i: 


(B-24) 


(B-25) 


(B-26) 


(B-27) 


(B-28) 


(B-29) 


(B-30) 


(B-31) 
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(B-32) 


dZ 


3L 


at 


0 




. VAR2 


(B-33) 


3Z 


4L 


= 0 -V > V > V 
at LIM L LIM 


(B-34) 


az 


5L ^ _±_ 

at Lg^ 


^AC t"4L - "5 l1 


S ^5L ^ ^ aAg ^2 ^5L aAg 
2 ] at n, 3t 


(B-35) 


The "A" matrix, 
given by; 



(i 


5, j = 1, . . .,5) for the upper pole is 


0 

0 

0 

^14 

0 

^21 

^22 

0 

^24 

0 

0 

^32 

0 

0 

0 

^41 

^42 

®43 

^44 

®45 

0 

0 

0 

®54 

^55 


(B-36) 


where 


a 


14 



(B-37) 
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(B-38) 


azi =1/K3 


"22 =-VS 


"24 = 


"32 


"41 


"42 K» p 


LE \ ‘'3 


"43 “ 


1 h \d 




^45 “ \c^^LE 
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(B-47) 




-ns fc^Lg 


^55 Lg \ 2 ^ ^AC 

° \ ^^3 


The following elements are zero: 


^11’ ^12’ ^13’ ^15’ ^23’ ^25’ ^31’ ^33’ ^34 ‘ ^35* ^51’ ^52’ ®53‘ 


The "A" matrix for the lower pole is given by: 


Elements through a^^ identical to the upper pole "A" matrix 

elements. The remaining two elements are: 


_ "4\c 

*54 Lg 




If V„ < -V.. .rw or V„ > the following upper pole elements are set to: 

U LIM U LIM o r iT 


**41 ■ °42 ■ “43 ■ ° 
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^44 ”^^AC 


(B-52) 


If V < “V or V > V the following lower pole elements are set to: 

L LiXM Li LiL^ 


^41 “ ®42 “ ° 


(B-53) 


"4.4 = -^\c 


LE 


(B-54) 


For steady state initialization: 


^lU ^2U ^IL ^2L ° 


(B-55) 


Zq__ “ Zq_ I R_ 
3U 3L o D 


(B-56) 


2y I T “ ^Ctt ~ 1 

4U 5U 4L 5L o 


(B-57) 
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APPENDIX C 


MATRIX INVERSION TECHNIQUE 


Given the matrix B 


f— 

— 

B , ! 


n-1 1 

u 

1 

n 

1 

V 1 

b 

n . 

nn 


t — 


(C-1) 


where is a (n-1) x (n-1) square matrix, is a 1 x (n-1) row, 

matrix, u„ is a (n-1) x 1 column matrix, b „ is a scalar. Assume the 
’ n ’ nn 

inverse to be of the form 





1 

C 1 

n-1 

1 



1 6-1 
1 n 


L 

1 _ 


(C-2) 


where is a (n-1) x (n-1) square matrix, q^ is a 1 x (n-1) row 

matrix, r^ is a (n-1) x 1 column matrix, 6 ^ is a scalar. Using the 
relationship BB ^ = I as in reference 2, it follows that: 


6 = b -V B^.u 

n nn n n-1 n 


(C-3) 


q 


n 



V B 
n 


-1 


n-1 


(C-4) 
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(C-5) 


r = -B 
n 


-1 

n-1 


u 

n 


6 


-1 

n 


C 


n-1 




u 

n 


q 


n 


(C-6) 


The 5x5 matrix Inverse was calculated by starting with n = 2 to get the 


2x2 inverse, with 

n = 3 to 

get the 3x3 inverse. 

etc. , 

until n = 5. 


If B is of 

the form 








^1 

0 

0 


0 




^21 

^22 

0 

^>24 

0 



B = 

0 

^^32 

^33 

0 

0 


(C-7) 


^1 

^2 

^3 

^4 

^5 




0 

0 

0 

^4 

^5 




— 




— 



the final 5x5 inverse, D, is given by 

the following equations 

. Equations (C-8) 

through (C-23) give the 4x4 

Inverse. 

Equations 

(C-24) 

through (C-33) 

give the 


5x5 Inverse extended from the 4x4 inverse. 


44 


^44 “ ^14 ^1 ^^24 ”^2 


(C-8) 


41 


= - d 


44 1 


(C-9) 


42 


= - d,, T„ 

44 2 


(C-10) 


22 



*43 


^44 ^43 
^^33 


(C-11) 


‘14 


*^44 ^14 


(C-12) 


I _ ^14 ^2 1 

24 44 i b^2 bjj ^22 


(C-13) 


I ^14 ^21 ^32 _ ^24 ^3 2 

34 44 \ b ^2 ^33 l »22 ^33 


(C-14) 


d 
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‘^41 ^3 


(C-15) 


1 = 21 

‘21 b^^ b22 41 ^4 


(C-16) 


^ ^21 ^32 _ 

'’ll ^^22 ^^33 ^ 


(C-17) 


^12 " ^42 ^3 


(C-18) 


‘'22 b^2 " ^^42 ^4 


(C-19) 
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(C-20) 


"32 ^ '^22^^33 - ‘" 42’^5 


“ *^43 "^3 


(C-21) 


^^23 “ ^43 3^4 


(C-22) 


^33 b33 " ^^43 "^5 


(C-23) 


55 b^3 - b^^ b^3 


(C-24) 


‘^51 “ *^55 ^54 *^41 


(C-25) 


*^52 “ ^55 ^54 ^*42 


(C-26) 


^^53 “ ‘^55 ^54 ^^43 


(C-27) 


^^54 *^55 '^54 *^44 


(C-28) 


24 



(C-29) 


‘15 


" ^55 ^^5 


25 


■ ^^55 \5 *^24 


35 


“ *^55 ^45 *^34 


"45 


" ‘^55 ^45 *^44 


(C-30) 


(C-31) 


(C-32) 


‘ij 


*^55 ^45 ^54 *^4j J ^ 


(C-33) 


where 



^42 ^21 ^43 ^21 ^32 

^11 ^22 ^11 ^^22 ^33 


(C-34) 


^42 _ ^43 ^32 

b22 b22 ^23 


(C-35) 


T 


3 



(C-36) 


T 


4 


^21 ^14 ^ ^2A 
^11 ^^22 ^22 


(C-37) 
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and 


^21 ^32 ^14 _ ^32 ^24 

^11 ^22 ^33 ^22 ^33 


(C-38) 
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TABLE V.- MBA PARAMETERS 


Maximum force output, N 

Bandwidth, m , rad/ sec 

O 

Bias current, 1^, A 

Nominal gap, g , 

N.m^ 

Actuator force constant, K, — ^ 

. A 

Gap coefficients, ^ and ^,/m 

Mass inductance measured at nominal gap, Lg^, H 

Leakage Inductance, L j p , H 

AC resistance, ^AC , ohms 

DC resistance, R^, ohms 


34.25 

628 

0.55522 

0.00762 

0.00161284 

127.795 

0.4347 

0.243 

238 


7.4 
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TABLE I.- 2.5 HZ SINE WAVE COMPARISON (BOX DENOTES MAXIMUM ERROR) 


t 

VOADAM 


210000000 

.117449790 

.148293576 

211000000 

.054503260 

.085357574 

212000000 

-.008456853 

.022400377 

213000000 

-.071415009 

-.040562478 

21A000000 

-.134355673 

-.103515452 

215000000 

-.197263310 

-.166443010 

216000000 

-.260122395 

-.229329620 

217000000 

-.322917413 

-.292159762 

218000000 

-.385632865 

-.354917929 

219000000 

-.448253272 

-.417588630 

220000000 

-.510763177 

-.480156398 

221000000 

-.573147152 

-.542605789 

222000000 

-.635389798 

-.604921389 

223000000 

-.697475750 

-.667087815 

224000000 

-.759389684 

-.729089723 

225000000 

-.821116316 

-.790911307 

226000000 

-.882640407 

-.852538807 

227000000 

-.943946772 

-.913955510 

228000000 

-1.005020275 

-.975146754 

229000000 

-1.065845839 

-1.036097435 


SF 3F 

LL, ^ 0 error, — = 0 

3 F 

error, ^ t 0 

a t 

.116824035 

.030043786 

-.000625755 

. 053876353 

.030854314 

-.000626907 

-.009084756 

.030857230 

-.000627904 

-.072043755 

.030852532 

-.000628745 

-.134985104 

.030840221 

-.000629431 

-.197893272 

.030820301 

-.000629961 

-.260752731 

.030792775 

-.000630336 

-.323547967 

.030757650 

-.000630555 

-.386263483 

.030714936 

-.000630618 

-.448883798 

.030664641 

-.000630526 

-.511393456 

.030606779 

-.000630279 

-.573777028 

.030541363 

-.000629876 

-.636019116 

.030468409 

-.000629318 

-.698104356 

.030387935 

-.000628605 

-.760017422 

.030299961 

-.000627738 

-.821743032 

.030204509 

-.000626716 

-.883265948 

.030101600 

-.000625540 

-.944570983 

.029991262 

-.000624211 

•1.005643003 

.029873520 

-.000622728 

•1.066466932 

.029748405 

-.000621093 




CO 

o 


TABLE II.- 10 HZ SINE WAVE COMPARISON (BOX DENOTES MAXIMUM ERROR) 


t 

VOADAM 

3F 

LL. ^ = 0 

3F 3F 

LL, — ^ 0 error, — = 0 

3t 3t 

dF 

error, — t 0 
* v.t 

,460000000 

.513067197 

.638874723 

.511264777 

.125787525 

-.001822420 

.461000000 

.256069556 

.382617966 

.254010125 

.126528409 


-.002079432 

.462000000 

-.001919442 

.124848986 

-.004248421 | 

1 

.126768427 

-.002328980 

.463000000 

-.259912634 

-.133413296 

-.262489983 

,126499338 


-.002577349 

.464000000 

-.516860116 

-.391147979 

-.519693729 

.125732137 


-.002813613 

.465000000 

-.771808873 

-.647336238 

-.774842931 

.124472635 


-.003034058 

.466000000 

-1.023668522 

-.900965369 

-1,026928999 

.122723153 


-.003240477 

,467000000 

-1.271522438 

-1 . 151032816 

-1.274955495 

.120489622 


-.003433056 

,468000000 

-1.514330397 

-1,396550151 

-1.517942083 

.117780246 


-.003611686 

. 46P000000 

-1.751152373 

-1.636546999 

-1,754928434 

.114605374 


-.003776061 

.470000000 

-1.981052266 

-1.870074899 

-1.964978028 

.110977368 


-.003925762 

.471000000 

-2.203121571 

-2.096211060 

-2.207181683 

.106910511 


-.004060312 

.472000000 

-2.416462944 

-2.314062036 

-2.420662156 

.102420908 


-.004179212 

,473000000 

-2.620293642 

-2.522767259 

-2.624575626 

.097526382 


-.004281985 

.474000000 

-2.813748851 

-2.721502461 

-2,818117040 

.092246390 


-.004368189 

.475000000 

-2,996084862 

-2.909482934 

-3.000522296 

.086601927 


-.004437435 

.476000000 

-3.166582075 

-3.085966639 

-3.171071474 

.080615435 


-.004489400 

.477000000 

-3.324567843 

-3,250257140 

-3,329091675 

.074310703 


-.004523832 

.478000000 

-3.469419126 

-3.4017C6355 

-3.473959679 

.067712771 


-.004540553 

,479000000 

-3.600564942 

-3.539717111 

-3.605104402 

.060847831 


-.004539460 





TABLE III.- TRIANGULAR WAVE COMPARISON (BOX DENOTES MAXIMUM ERROR) 


t 

VOADAM 

TT 

LL, — = 0 
3t 

3F 9F 

LL. — ^ 0 error. — = 0 

3F 

error, -r— ^ 0 
0 t 

5A1000000 

.499271067 

.498859295 

.499337367 

.000411772 


. 000066300 

542000000 

.499841323 

.499663338 

.500043179 

.000177985 


.000201856 

543000000 

. A99947785 

.500023750 

.500181731 

.000075964 


.000233946 

544000000 

.499615176 

. 499900124 

.49982804 8 

.000292949 


.000212873 

545000000 

.490926654 

.499380925 

.499095039 

.000454271 


. 000168384 

546000000 

.497978622 

.498539647 

.498096720 

.000561025 


.000118098 

547000000 

.496860664 

.497482348 

.496931430 

.000621683 


.000070765 

548000000 

.495645670 

.496293016 

.495676123 

.000647145 


.000030253 

549000000 

.494388219 

.496036759 

.494386464 

.000648540 

- 

.000001755 

550000000 

.493124903 

.493760065 

.493099624 

.000635162 

- 

•.000025279 

551000000 

.491879290 

.492493394 

.491838045 

.000614105 

- 

•.000041245 

552000000 

.490664179 

.491254554 

.490613202 

.000590375 

- 

-.000050977 

553000000 

.489484820 

.490052042 

.489428937 

.000567222 

- 

-.000055883 

554000000 

.488341439 

.488887967 

.468284159 

.000546528 

- 

-.000057281 

555000000 

.487231212 

.487760394 

.487174910 

.000529182 

-.000056302 

556000000 

.486149720 

.486665122 

.486095856 

.000515402 

-.000053865 

557000000 

.485091961 

.485596948 

.485041292 

.000504986 

-.000050669 

558000000 

.484053002 

.484550510 

.484005779 

.000497508 

-.000047223 

559000000 

.483028362 

.483520807 

.482984495 

.000492445 

- 

-.000043867 

560000000 

.482014203 

.482503472 

.481973390 

.000489269 

. 

-.000040813 





TABLE IV.- STEP COMPARISON (BOX DENOTES MAXIMUM ERROR) 


t 

VOADAM 

TT 

LL, — s 0 

at 

error, = 0 

0 t 

005000000 

0.000000000 

0.000000000 

0,000000000 

006000000 

,162252170 

.147347430 

-.014904740 

007000000 

.493552465 

.480234865 

-.013317600 

000000000 

.84416691? 

.637443361 

-.006725552 

000000000 

1.146302109 

1.146273165 

-.000028943 

010000000 

1.376957376 

1.361929791 

.004972415 

011000000 

1.536057579 

1.543985964 

.007928385 

012000000 

1.633690697 

1 .642993667 

.009102990 

013000000 

1.68416061P 

1.693122399 

.008961780 

OIAOOOOOO 

1.700371363 

1.706349390 

.007978028 

015000000 

1.694159157 

1.700710234 

.006551077 

016000000 

1.674726022 

1.679707480 

,004981459 

017000000 

1,648861594 

1.652335073 

.003473479 

018000000 

1.621246513 

1.623396321 

.002149808 

019000000 

1.594659358 

1.595929236 

.001069880 

020000000 

1.571368877 

1.571636766 

,000247909 

021000000 

1.551602062 

1.551270672 

-.000331410 

022000000 

1.535648199 

1.534947799 

-.000700399 

023000000 

1.52329AQ06 

1.522395924 

-.000898982 

02<»000000 

1.514102440 

1.513134309 

-.000968131 
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Figure 1.- Magnetic Bearing Assembly Model 
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Figure 2.- 


2.5 Hz sine wave comparison (+ symbol 
denotes VOADAM solution) 









Figure 3.- 10 Hz sine wave comparison (+ symbol 
denotes VOADAM solution) 
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